function PD = PD(p,gH,gL,r)


PD = p*Gordon(gH,r) + (1-p)*Gordon(gL,r);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


function P = Gordon(g,r)

if g>r
    error('gordon growth formula invalid')
end

P = (1+g)/(r-g);





